Magnetic couplings, optical spectra, and spin-orbit exciton in 5d electron Mott 

insulator Sr 2 Ir04 
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Based on the microscopic model including spin-orbit coupling, on-site Coulomb and Hund's in- 
teractions, as well as crystal field effects, we have investigated magnetic and optical properties 
of Sr2lrC>4. Taking into account all intermediate state multiplets generated by virtual hoppings 
of electrons, we calculated the isotropic, pseudodipolar, and Dzyaloshinsky-Moriya coupling con- 
stants, which describe the experiment quite well. The optical conductivity o~{uj) evaluated by the 
exact diagonalization method shows two peaks at ~ 0.5 and ~ 1.0 eV in agreement with experiment. 
The two peak structure of cr(w) arises from the unusual Fano-type overlap between electron-hole 
continuum of the J e // = 1/2 band and the intrasite spin-orbit exciton observed recently in Sr2lr04. 

PACS numbers: 75.30.Et,71.70.Ej,78.20.Bh 



Mott physics is one of the most fundamental phenom- 
ena in condensed matter physics giving rise to diverse and 
fascinating collective behavior of correlated electrons [l| . 
In Mott insulators, strong on-site Coulomb repulsion (U) 
splits the half-filled band into the lower Hubbard band 
(LHB) accommodating the spin and orbital degrees of 
freedom of electrons, and the empty upper Hubbard 
bands (UHB). The electrons hop into adjacent sites only 
virtually, overcoming thereby the Mott-Hubbard gaps. 
Typically, the antiferromagnetic (AFM) ground state is 
realized due to the kinetic energy gain of the virtual ex- 
change processes. The optical conductivity shows peak 
structures at the absorption edges corresponding to the 
transitions between LHB and UHBs. 

Recently, a new class of 5d Mott insulators such as 
Sr 2 Ir04 and Na2lr03 0, @ has been discovered, 

where the strong spin-orbit (SO) coupling is crucial for 
stabilizing the insulating state. The Coulomb interaction 
between 5d-electrons U ~ 2 eV is much smaller than that 
in conventional 3d electron Mott insulators. On the other 
hand, the SO coupling A ~ 0.4 eV in iridates is much 
larger and splits 2 T l£( states of 5d 5 -shell into the half- 
filled J e ff = 1/2 and fully-occupied J e // = 3/2 states. 
Then the narrow J e ff = 1/2 band tends to undergo Mott 
transition even at relatively small U 0, 0, 0-E1 • 

The issue in iridates is to which extent the physical 
properties and model descriptions of J e // = 1/2 sys- 
tems are (dis) similar to those of conventional 3d Mott 
insulators. The best studied member of 5c? Mott sys- 
tems is Sr2lr04 perovskite. Its in-plane canted AFM 
ground state [H , magnon spectra [S5| , and finite temper- 
ature spin dynamics [11| closely resemble those of parent 
high-T c cuprates, in accord with the theoretical predic- 
tions [12j based on Mott-Hubbard picture. In addition to 
magnons, resonant inelastic x-ray scattering (RIXS) ex- 
periments 13, si| have observed also the higher energy 
broad peak at ~ 0.5 — 0.8 eV. Based on the theoreti- 
cal expectations 12, 3 IH that the SO split t2 g man- 



ifold should have a magnetically active mode at ~ |A, 
this peak has been attributed to the transition between 
J e ff = 1/2 and 3/2 states and termed "SO exciton" [sBj ]. 

Concerning the charge excitation spectra, the optical 
conductivity cr(ui) of Sr2ir04 shows two peaks at ~ 0.5 
and ~ 1.0 eV [2j that are preserved up to high tempera- 
tures (lq |. The first peak was assigned to the transition 
from occupied J eff = 1/2 LHB to J eff = 1/2 UHB, 
while the peak at ~ 1.0 eV to that from J e ff — 3/2 to 
unoccupied J e ff = 1/2 [J . This interpretation, however, 
is based on a picture of single-electron density of states. 
In fact, there has been no theoretical calculations of oiiS) 
in iridates taking into account the many-electron multi- 
plet structure of excited states, which is known to be es- 
sential for the interpretation of optical data. Moreover, 
there are inherent relations between the optical absorp- 
tion peaks/intensities and the strength of magnetic cou- 
plings (both are determined by the same virtual hoppings 
and excited multiplets), which enable one to extract the 
physical parameters, such as U and Hund's coupling Jh, 
from a combined analysis of the magnetic and optical 
data (l7rfl9| . The aim of the present Letter is to extend 
this fruitful approach to the 5d electron Mott insulators. 

We have calculated the magnetic couplings, optical 
conductivity, and RIXS spectra in Sr 2 Ir04 by exact di- 
agonalization (ED) of a microscopic model on small clus- 
ters, fully incorporating the multiplet structure of Ir ions, 
SO coupling, tetragonal distortion, and octahedral rota- 
tions. The magnetic co uplings obtained are consistent 
with the available data @,[lH[S5j]. Calculated RIXS spec- 
tra reproduce the SO exciton mode. More interestingly, 
we found that, unlike the case of 3d oxides [l7l4l9j . the 
observed peaks in <t(uj) of Sr2lr04 cannot be directly de- 
termined from the multiplet energies and intensities. In- 
stead, the strong mixing between intersite optical excita- 
tions (electron-hole continuum) and intraionic transition 
between J e ff = 1/2 and 3/2 states (SO exciton) is found 
to be essential feature of Sr2ir04. This implies that nei- 
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ther purely atomic nor simple band picture is sufficient 
to describe the small charge-gap iridium oxides. 

Model- To describe the electronic structure of Ir ions, 
we adopted the following Hamiltonian: 
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where r and isospin a [12| refer to the lowest three 
Kramers doublets [the eigenstates in the presence of the 
tetragonal crystal field splitting A xy and SO coupling, see 
Fig.QJa)]. The other terms describe the on-site Coulomb 
and Hund's interactions, where u, and a are the orbital 
and spin indices, respectively. As usual, we parametrize 
the interaction matrix as = U, V = U — 2J#, 
and J„ v =J'^ = JhM- 

Fig. [ljb,c) shows all possible multiplets of c? 5 and 
d 4 electronic configurations. For d 6 , all three lowest 
Kramers doublets t 1i2j3 are fully occupied. We note 
that the wave-functions of these states are often taken 
in the limit of lODq — > oo, i.e., admixture of e g or- 
bitals (oc X/lODq) in the ground state via SO coupling 
is ignored. Since the SO coupling in iridates is strong, 
we fully include the e g orbital admixture in the lowest 
Kramers doublets. 

Because an electron hops between nearest-neighbor 
(NN) Ir sites via an oxygen [see Fig. Eel)], and the 
charge transfer energy A between Ir-5d and 0-2p is large 
(A ~ 3.3 eV [IH), we assumed the following effective NN 
hopping Hamiltonian: 
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where the effective hopping integral between tct and t' a' 
states is t iT . 5 . 7T / 5 / = V* t " a *<' t iT' „'■,?<, calculated 

the pd-hopping matrix ij T5 . ;p(T between 5d Ir and 2p O 
orbitals in terms of two parameters t p d a and t p dn 23] . 

Magnetic interactions - We considered the Iri-Ir2 pair 
described by the following Hamiltonian: 



H — Hi + Hi 
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Here Hi and Hi are the ionic Hamiltonians given by 
Eq. (J]), and Hu is the hopping term of Eq. ([2]). Us- 
ing the ED method, we solved this Hamiltonian nu- 
merically, taking into account all possible multiplets al- 
lowed for the d 5 — d 5 and d 4 — d 6 configurations, and 



TABLE I: Physical parameters in units of eV. 
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FIG. 1: (color online) (a) Energy levels of d electron in 
the presence of the cubic crystal field lODq, the tetragonal 
crystal field A. xy , and the SO interaction A [2j]. The lowest 
multiplet levels of (b) d 5 and (c) d 4 configurations including 
SO A and the Hund's coupling Jh- (d) Two-site and (e) 
four-site clusters in the xj/-plane, which were employed in the 
evaluation of the magnetic couplings and conductivity u(u). 



obtained eigenvalues E n and eigenstates \ip n ) of the Ir- 
Ir cluster. The Hamiltonian can be then expressed as 
H = J2 n -^nlV'nXV'nl- Without the hopping term Hu, the 
Ir-Ir pair would have four degenerate states denoted by 
|^>°) (i — 1-4). The magnetic interactions are generated 
by virtual hoppings among them. The resulting effective 
magnetic Hamiltonian can be obtained by applying the 
projection operator 7 , i/2 = X]i=i l^iX^i I on t° H 24|: 



Vi/iHTi/i - -Tv (Vi/iHVi/i) =Si-J-S 2 , (4) 

where S is isospin one-half, and J is a 3 x 3 tensor. 

Consistent with symmetry considerations as well as 
with Ref. [lH, we found that Eq. (Q} comprises four dis- 
tinct nonvanishing terms, 

JSy S 2 + 5J z Si >z S 2 , z + SJ xy (Si-rn)(Si- fi 2 ) + D • Si x S 2 , 

(5) 

corresponding to isotropic Heisenberg (J), symmetric 
{5J Z , SJ xy ) and antisymmetric (D) anisotropic couplings 
between NN Ir-moments. We note that only the z- 
component Dzyaloshinsky-Moriya (DM) vector D sur- 
vives because of the mirror symmetry with respect to 
the xy-plane, i.e., D = (0,0,1?). 

Figure[5]presents the magnetic coupling constants, cal- 
culated by using the parameters provided in Table HI as 
a function of Ir-O-Ir bonding angle <j >. 10-Dg and A 



values are adopted from Refs. [13|, 122J, [25J and A from 
Refs. [13, EH- We have determined the optimal values 
of U, Jh, t p d such that they describe both the magnetic 
and optical data properly [27(. We set tetragonal split- 
ting A xy — 0.15 eV; at this value, |D/J| ratio becomes 
~ 0.34 and yields the spin canting angle of ~ 9.3° at 
<j> = 158° (as in Sr 2 Ir04). As shown in the inset of Fig. [21 
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FIG. 2: (color online) Magnetic coupling constants between 
NN Ir ions (r*i2 1 with respect to the bonding angle (<j>). J xx , 
represent diagonal parts of the superexchange 
; = J + 8J xyi Jyy = J, and J zz = J + SJ Z are 



J yy i SUld Jzz 

tensor J . J x 



deduced from Eq. ((5j). D is a z-component of the DM vector 
D. Arrows denote the bonding angle in SrfrCU. Inset shows 
\D/J\ ratio as a function of the tetragonal splitting A xy . 
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FIG. 3: (color online) (a) Schematic diagrams of possible 
multiplets included in subsets (Hi) of the Hilbert space. D, 
Q, S, and T refer to multiplets of d 5 and d 4 , as labeled in 
Fig-Hb) and (c), and A denotes a nondegenerate d 6 state. H5 
includes also diagrams (not shown) with P,P',S' multiplets 
of the d 4 state [Fig. [2(c)] instead of T. (b) Two examples 
of optically active transitions that contribute to <r(w) at low 
energies. The final states H3 and H5 have an overlap with on- 
site local SO exciton Q (i.e., with H2 sector), due to intersite 
hoppings between J e ff = 1/2 and J e ff = 3/2 states, resulting 
in two peak structure of o(oj). 



when A xy is too small, \D/J\ ratio becomes too large to 
describe the observed spin canting. Approximately, we 
find that D oc sin (2</>). The spin canting angle of ~ 9.3° 
is close to those found in previous studies [l2|, HH . Note- 
worthy is that the Ising coupling SJ Z (= J zz — J yy ) in 
Eq. ([5]) is enhanced when decreasing bonding angle <fi. 
We also calculated a pseudodipolar coupling 5 J xy for dif- 
ferent Jh values and found that SJ xy becomes zero at 
Jh = 0. Both of these observations agree with the ana- 
lytical results (l2| . 

Summarizing our results for magnetic interactions in 
Sr2lrC>4, we obtained the following values at 4> = 158°: 
J ~ 76.8, D ~ -26.2, 6J Z ~ 5.9, and SJ xy ~ 8.6 meV. 
With these coupling constants, the interactions of Eq. ([5]) 
lead to the canted AFM state with Ir-moments lying in 
the zy-plane, as observed in Sr 2 Ir04 \$. The calculated 
isotropic coupling J is in close agreement with the ex- 
perimental value of J ~ 60 meV |S5| . Unusually large 
anisotropic couplings D,5J z ,5J xy found here are the di- 
rect fingerprints of strong SO interaction. 

Optical conductivity and RIXS spectra.- In 3d Mott 
insulators, the hopping between NN sites plays a dom- 
inant role in o~(uj) that shows peak structures near the 
ionic multiplet states — 17-fl9j|. and 2-site clus- 
ter may capture essential features of <j(lo). In contrast, 
5d Mott insulators have weaker Coulomb repulsion and 
thus a duality of atomic and band nature of correlated 
electrons is more pronounced S5[. In order to capture 
the derealization of optically excited electron-hole (e-h) 
pairs, we consider here a 2 x 2 cluster shown in Fig. [TJe). 

We considered all possible multiplets within the d 5 -d 5 - 
d 5 -d 5 and d 4 -d 6 -d 5 -d 5 charge configurations. In order to 



clarify the origin of optical peaks, it is useful to classify 
the Hilbert space into 6 subspaces [see Fig. [3fa)]: Hi of 
four d 5 doublets (D), H 2 of one or more quartets (Q) 
among four d 5 configurations, H 3 of two D's of d 5 and 
NN e-h pairs d 4 (S)-d 6 (A), H 4 of two D's of d 5 and next 
nearest-neighbor (NNN) e-h pairs d 4 (S)-d 6 (A), H 5 of two 
D's of d 5 and (NN & NNN) e-h pairs d 4 (T, P, P' , S')- 
d 6 (A), and the remaining states (He,...) that involve, 
e.g., simultaneous intersite e-h transitions and local SO 
exciton Q. All the above configurations couple to each 
other via the hopping Hamiltonian of Eq. © . 

We have solved the Hamiltonian matrix with the ED 
method, and obtained <j(oS) from the following relation: 

a(uj) = 7TV Pn\(tp m \Jc\lpn)\ 2 8(uJ + E n —E m ), 

(6) 

where v is volume per Ir-site, p n is the probability density 
of eigenstate \ip n ), and J c is the current operator. We set 
B^ 1 = ksT = 30 meV to avoid the finite size effect. A 
function S(u) is treated with broadening of 0.05 eV. 

Figure EJa) shows the result for a(oj) calculated by us- 
ing the parameters from Table HI Two peaks are revealed 
clearly at ~ 0.5 and ~ 1.0 eV, in good agreement with 
experiment. Estimated value of <r(uS) at ~ 0.5 eV is also 
consistent with the experimental data. For lu > 0.6 eV, 
calculated <r(ui) is smaller than that observed. The pos- 
sible reasons for this discrepancy are the contributions 
from the pd-charge transfer peak at ~ 3 eV and from 
two or more e-h pair excitations not included here. 

Shown in Figure IDJb) is the RIXS spectra, calculated 
by employing the inelastic x-ray scattering operator of 
Ref. [15| instead of J c in Eq. ©. At the Z*3-edge, intense 
magnon (below 0.25 eV) and SO exciton (0.5 — 1.0 eV) 
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FIG. 4: (color online) (a) Optical conductivity calculated 
(solid line) and experimental data [2| (dashed line). Red ver- 
tical sticks show relative strengths and positions of optical 
transitions without broadening. Inset: the result for ct(lj) 
when the engenstates \ip n ) and energies E n in Eq. (|6]) are ap- 
proximated by purely local ionic multiplets. (b) Calculated 
Z/3-edge (solid line) and 10 x Z/2-edge (dotted line) RIXS spec- 
tra at q = (n,n). Dashed line: experimental data |S5| . (c) 
Projected excitation spectra: Ai represents magnon band, A2 
contains one or more SO excitons, A3 and A4 represent NN 
and more distantly separated e-h excitations derived from the 
J e ff = 1/2 states, respectively, and A 5 shows the e-h contin- 
uum of J e ff = 3/2 states. A6 refers to other excitations, e.g., 
simultaneous transitions in e-h and SO exciton channels. 



bands are obtained, while the intensity at the L2-edgeis 
nearly vanishing, in agreement with experiments [T3L IS5jj . 

Of particular interest is the origin of two peaks in u(uj) 
and their relation to SO exciton. In 3d Mott insulators, 
each peak would correspond to the specific multiplet of 
eP" 1 — d n+1 with th e sp ectral weight proportional to 
I (d n - 1 (f 1+1 \J c \d n d n )\ 2 [lM9| . Such a simple ionic multi- 
plet picture, however, is incomplete in iridates. Indeed, 
as shown in inset of Fig. @|a), a(u>) based on this pic- 
ture shows that the first peak [\DDDD) -> \ASDD) 



transition in Fig. [3[b)] is much stronger than those cor- 
responding to higher energy transitions involving T, P- 
multiplets [29|, in contrast to the observation that the 
spectral weight of ~ 1.0 eV peak is even larger than that 
of -0.5 eV peak [Fig. Ufa)]. 

To understand how the e-h derealization effects lead 
to strong deviations from a simple ionic picture, it is 
instructive to analyze the underlying excitation spectra 
in more detail. To this end, we evaluated the projected 
excitation spectrum (PES) of the 2x2 cluster into the 
subspaces introduced above: 

A *M=E E \(^nH\ 2 S(u-E n ), (7) 



where \m) represents the orthonormal basis of the sub- 
space Hi. We can identify the PES in Fig.^c) as follows: 
Ai represents magnon sector, A2 is related to one or more 
intra-site D — Q transitions (SO excitons), A3 and A4 are 
the e-h continuum of J e // = 1/2 states, and A 5 describes 
the e-h continuum of J e ff = 3/2 bands. 

We notice intriguing features in Fig. |4jc): (i) Ai and 
A2 show the peaks in the range of — 0.25 and 0.5 — 1.0 
eV, respectively. These peaks are manifested in the RIXS 
spectra as the magnon and SO exciton bands. A2 shows 
also the high energy peak at — 1.5 eV, which however 
does not directly couple to the RIXS process since it 
contains two SO excitons residing on different sites, (ii) 
A3 and A4 (J e ff = 1/2 e-h continuum) are located in the 
wide range of 0.5 — 1.5 eV. This fact reveals that not only 
the lower but also the higher peak of <r(u>) is attributed to 
the e-h continuum of J e f / = 1/2 states, contrary to pre- 
vious interpretation based on a simple band picture 0. 
(iii) A3 (NN e-h contribution) is depleted in the vicin- 
ity of the SO exciton A2. The mixing of these modes 
is natural in view of the nonzero overlap between quasi- 
degenerate \ASDD) € H 3 and \DQDD) £ H 2 states [see 
Fig. [3ja,b)] by virtue of intersite hoppings. 

The PES thus shows that the optically active e-h con- 
tinuum of Jef f = 1/2 band and the optically forbidden 
SO exciton [S5J are located in the same energy range of 
0.5 — 1.5 eV, and there is considerable mutual interaction 
between these excitations due to the interband hopping 
between J e ff = 1/2 and J e // =3/2 states. Two peak 
structure in a(u>) results from this unusual mixing among 
different excitations, which is reminiscent of the charac- 
teristic behavior of the Fano resonance. 

To conclude, we presented a unified description of mag- 
netic couplings, optical conductivity, and RIXS spectra 
in Sr2lr04 within the model including strong SO cou- 
pling and ionic multiplet effects. The results obtained are 
consistent with the available experimental data. Since 
the SO splitting and the Mott-Hubbard charge gap are 
of similar scale, an unusual coupling between two dif- 
ferent types of excitations - the optically inactive SO 
exciton and the optically active e-h continuum - is in- 
duced due to the interband mixing of J e ff = 1/2 and 
J e ff = 3/2 states. Although this effect is not very es- 
sential for NN magnetic interactions, it plays a crucial 
role in determining the shape of optical spectra: the e-h 
continuum is suppressed in the vicinity of the SO exci- 
ton inherent to iridates. The unusual Fano- type cou- 
pling between the rich ionic multiplet excitations and 
electronic continuum seems to be a characteristic fea- 
ture of Jeff = 1/2 Mott insulators. This phenomenon 
is rooted to the atomic/band duality of 5d electrons and 
should therefore be generic to iridates with unusual hier- 
archy of energy scales. 

We thank G. Jackeli for useful discussions. This work 
was supported by the NRF (No.2009-0079947). 
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SUPPLEMENTAL MATERIAL 



The gap value and peak positions in optical spectra 
are most sensitive to the interaction parameters U and 
Jh- At the same time, there is large uncertainly con- 
cerning their values in iridates: current estimates for U 
(~ 2.2 eV [HIH, 3 eVjsl) and for J H (~ 0.2 eV [§3], 
0.3 eV [iH, 0.6 eV [si[sj) vary broadly. We therefore 
determine U and Jh values such that they give a reason- 
able agreement with the observed optical gap and peak 
positions, and, at the same time, we cross-check them by 
calculating the magnetic exchange constants. 

We found that U ~ 2.2 eV and J H ^ 0.2 - 0.3 eV sug- 



gested by Refs. [Sl|,lS2| result in too large optical gap (see 
a representative curve in Fig. IS1|) and too small magnetic 
exchange constant J. Far better agreement with experi- 
ment is obtained when we slightly descrease U and sub- 
stantially increase Jh (above 0.4 eV). Solid line in Fig.lSlI 
shows the result for U = 1.86 eV and Jh = 0.5 eV. As 
discussed in the main text, this parameter set provides 
also the magnetic interactions and spin canting angles 
that are in fair agreement with experiment. 

U and Jh can still be varied within the ranges 1.6-2.0 
eV and 0.4-0.5 eV, respectively, with acceptable fit results 
concerning the optical data. For instance, we show in 
Fig.EDthe result for U = 1.60 cV and J H = 0.42 eV (esti- 
mated as Jh = 3B+C using the Racah parameters B and 
C for an octahedrally coordinated Ir 4+ impurity (S6|). 
While these parameters provide optical spectra similar 
to that we found for U = 1.86 eV and Jh = 0.5 eV, the 
latter choice gives a magnetic coupling J which is closer 
to the observed value. 



* bimin@post ech.ac.kr| 
[SI] C. Martins, M. Aichhorn, L. Vaugier, and S. Biermann, 

Phys. Rev. Lett. 107, 266404 (2011). 
[S2] R. Arita, J. Kunes, A.V. Kozhevnikov, A.G. Eguiluz, and 

M. Imada, Phys. Rev. Lett. 108, 086403 (2012). 
[S3] R. Comin, G. Levy, B. Ludbrook, Z.-H. Zhu, C.N. Veen- 

stra, J. A. Rosen, Yogesh Singh, P. Gegenwart, D. Strieker, 

J.N. Hancock, D. van der Marel, I.S. Elfimov, A. Damas- 

celli, larXiv:1204.447Tl 
[S4] D. van der Marel and G.A. Sawatzky, Phys. Rev. B 37, 

10674 (1988). 

[S5] J. Kim, D. Casa, M.H. Upton, T. Gog, Y.-J. Kim, J.F. 
Mitchell, M. van Veenendaal, M. Daghofer, J. van den 
Brink, G. Khaliullin, and B.J. Kim, Phys. Rev. Lett. 108, 
177003 (2012). 

[S6] B. Andlauer, J. Schneider, and W. Tolksdorf, Phys. Stat. 
Sol. B 73, 533 (1976). 



t 1 r 




0.0 0.5 1.0 1.5 2.0 

co(eV) 

FIG. SI: (color online) Optical conductivity calculated with 
three different sets of U and Jh interaction parameters: U = 
1.86 eV and J H = 0.50 eV (solid line), U = 1.60 eV and J H = 
0.42 eV (dotted line), U = 2.25 eV and J H = 0.2 eV (dashed- 
dotted line). Green dashed line represents the experimental 
data. The Heisenberg exchange couplings J calculated with 
these parameter sets are J = 76.8 meV, 83.0, and 37.8 meV, 
respectively; the experimental value is J ~ 60 meV [S5l ]. 



